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ABSTRACT 

This  report  identifies  potential  errors  in  computing  high  statistical  reliability  for  a  required  com¬ 
ponent  fatigue  life.  The  reliability  values  were  deternjined  from  application  of  a  joint  probability  density 
(JPD)  analysis  used  in  an  American  Helicopter  Society  round  robin  safe  life  problem. 

In  the  analysis  normal  probability  density  functions  (PDFs)  were  assumed  for  both  the  material 
strength  and  the  spectrum  load  values.  The  PDF  model  parameters  were  varied  and  the  PDFs  were 
slightly  modified  (contaminated)  in  order  to  examine  the  sensitivity  in  computing  high  statistical  reliability 
when  uncertainties  exist  in  assuming  the  PDF.  Lower  tails  of  the  PDFs  were  also  modified  by  trunca¬ 
tion,  independent  of  the  model  contamination,  in  order  to  determine  the  relative  influence  on  reliability 
from  tail  modifications  as  compared  with  the  parameter  uncertainties  and  contamination.  The  stability 
of  statistical  estimates  of  the  extreme  tail  quantiles  and  their  corresponding  probabilities  as  a  function  of 
sample  size  were  examined  for  a  generic  distribution. 

Assuming  a  PDF  to  represent  load  or  material  strength  is  a  substantially  more  critical  issue  than 
accurate  representation  of  the  extreme  lower  tail  of  the  PDF  when  computing  high  reliability.  Sampling 
trials  for  extreme  tail  quantiles  and  reliabilities  indicate  that  unstable  values  can  result  from  sample  sizes 
of  100. 

The  primary  conclusion  from  these  analytic  results  is  that  the  computation  of  a  high  statistical 
reliability  may  have  little  or  no  association  with  actual  engineering  high  reliability. 
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INTRODUCTION 


The  use  of  a  quantitative  high  reliability  requirement  for  a  helicopter  component  fatigue  design 
has  received  considerable  attention  recently.  The  U.  S.  Army  established  a  requirement  of  .999999 
(“six  nines”)  reliability  for  dynamic  components  in  its  most  recent  helicopter  development1.  Sub¬ 
sequently,  the  American  Helicopter  Society  Subcommittee  on  Fatigue  and  Damage  Tolerance,  con¬ 
ducted  a  round  robin  study  of  high  reliability  fatigue  methodology  applied  to  a  simple  structural 
element2.  A  review  of  the  round  robin3  noted  difficulties  in  the  reliability  analysis.  Each  participant 
used  a  different  fatigue  curve  and  fatigue  limit  variability  which  resulted  in  significantly  different 
fatigue  lives,  for  the  six  nines  reliability  requirement.  A  recent  fatigue  analysis  by  a  helicopter 
manufacturer4  found  that  “  reliability  is  very  sensitive  to  changes  in  the  population  mean  strength 
and  scatter”.  In  addition  Reference  4  notes  “the  conclusions  of  this  study  are  not  fully  applica¬ 
ble  to  actual  fleet  management  due  to  the  presence  of  statistically  indeterminant  variables  such  as 
degraded  or  non-conforming  components”. 

The  present  authors  in  previous  study,  have  investigated  the  sensitivity  of  high  reliability  com¬ 
putations  from  a  stress-strength  model5  to  uncertainties  in  the  identification  of  the  probability 
density  functions(PDFs)  in  the  model.  The  uncertainties  are  associated  with  the  selection  of  com¬ 
peting  parametric  forms(  e.g,  normal,  log-normal,  Weibull,  etc.)  or  with  the  undetected  presence 
of  contaminated  populations.  Contaminated  distributions  could  be  bimodal,  caused  by  degraded 
or  non-conforming  components,  or  could  be  the  result  of  by  unexpected  loading  anomalies.  The 
results  from  Reference  5  showed  that  high  reliability  estimates  can  vary  substantially  even  for  “al¬ 
most  undetectable”  differences  in  the  assumed  stress  and  strength  PDFs.  The  authors  have  also 
investigated  the  sensitivity  of  safe  life  fatigue  reliability  of  a  simple  structural  element  loaded  by 
a  simplified  spectrum  to  a  variety  of  uncertainties6,  demonstating  that  a  small  amount  of  uncer¬ 
tainty  in  the  parameters  of  the  load  or  strength  PDF  resulted  in  a  substantial  reduction  in  the  high 
statistical  reliability  values  for  a  specified  lifetime  of  the  component. 

The  round  robin  review,  Reference  3,  also  expressed  a  concern  for  the  effect  of  inaccuracies  or 
truncations  of  the  tails  of  the  distributions.  An  investigation  of  the  truncation  of  known  normal 
PDF  was  proposed  in  order  to  determine  an  “acceptable  degree  of  truncation”  in  computing  high 
statistical  reliability.  Apparently,  this  determination  would  be  expected  to  indicate  the  portion  of 
the  tail  region  in  which  an  accurate  representation  of  the  PDF  is  not  required. 


In  this  report  the  AIIS  round  robin  fatigue  problem  and  its  methodology,  Reference  2,  will 
be  used  to  investigate  the  cited  issues  by  considering:  a)  The  effect  of  small  changes  in  the  PDF 
parameters  on  the  reliability-life  relationship,  b)  The  influence  upon  reliability  of  the  consideration 
of  PDFs  which  are  contaminated,  using  the  methods  of  Reference  5,  by  bimodal  effects,  c)  The 
“true”  reliability  associated  with  fatigue  lives  which  have  been  obtained  by  satisfying  an  “apparent” 
six  nines  reliability  based  on  normal  PDFs  which  have  been  truncated  in  the  extreme  tail  region, 
d)  The  relative  influence  on  high  statistical  reliability  of  parameter  uncertainties,  contamination 
and  truncations  of  the  PDF.  The  consideration  of  issues  involving  the  extreme  tails  of  the  PDFs 
requires  an  accurate  measure  of  the  truncation  point  locations,  which  is  difficult  to  achieve,  since 
sufficient  amounts  of  data  is  usually  not  available.  In  practice,  truncation  point  locations  would  be 
estimated  from  small  data  sets  of  load  or  strength  measurements.  The  stability  of  the  statistical 


estimates  of  the  extreme  tail  quantile  and  probability  values  will  be  investigated  based  on  sampling 
simulations  of  a  generic  normal  PDF. 

These  results  will  be  assessed  to  indicate  the  potential  role  of  the  PDF  tail  truncation  analysis 
in  providing  conclusive  information  on  the  acceptability  of  PDF  modeling  and  whether  a  quantified 
.96  reliability  provides  a  meaningful  measure  which  correlates  with  levels  of  structural  integrity. 


FATIGUE  LIFE  COMPUTATIONS 


The  following  standardized  fatigue  life  computation  procedures  were  obtained  from  a  round 
robin  study  conducted  by  the  AIIS,  Reference  2.  The  form  of  the  S-N  curve  is, 

N  =  C{S'  -  Se)d,  (1) 

where  N  =  number  cycles  to  failure;  Se  =  fatigue  strength  for  very  large  N  values,  for  mini¬ 
mum  stress  equal  to  zero;  S'  =  effective  maximum  cyclic  stress,  for  minimum  stress  equal  to 
zero, equivalent  to  spectrum  stresses;  C  and  D  are  parameters  from  regression  least  squares  analysis. 


In  order  to  apply  the  S-N  curve  in  Equation  1  using  the  actual  operating  load  spectrum,  the 
following  relation  for  S'  is  required: 


5* 


Q  ■  Su-  SL _ 

Su  -  a  ■  Sm  +  a  ■  Sl/2 


(2) 


This  equation  represents  a  form  of  the  Goodman  correction  factor  used  in  Reference  2,  which 
converts  a  defined  spectrum  mean  stress  and  stress  range  to  an  equivalent  stress  range  which  causes 
equal  fatigue  damage  from  zero  to  specified  S'  value.  Su  represents  the  ultimate  strength  of  the 
material.  Sm  and  Se  represent  the  mean  and  range  respectively  of  the  nominal  stress  from  a  rainflow 
count  obtained  in  Reference  2,  of  the  standardized  Felix  28  spectrum  as  tabulated  in  Table  1.  The 
a  value  is  a  scaling  parameter  for  the  spectrum  load  values  Si  and  Sm  representing  the  effective 
load  scaling,  over  the  lifetime  of  a  component.  This  parameter  can  provide  changes  in  the  baseline 
spectrum  load  in  order  to  account  for  differences  in  usage,  pilot  technique,  weather,  weight,  etc. 


Let  the  fatigue  life  Np  represent  the  number  of  passes  prior  to  the  component  failure.  Then 
from  Miner’s  Rule,  Reference  6, 


where, 


NP  =  1  /DF 

DF=Y^. 

h  mk) 


(3) 

(4) 


The  n(k)s  represent  the  number  of  cycles  for  a  specific  k  value  and  NK  represents  the  total  number 
of  spectrum  load  values  from  Table  1.  The  N(k)s  are  the  results  from  Equation  1. 
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RELIABILITY  AS  A  FUNCTION  OF  LIFE 


The  following  procedures  suggested  in  Reference  1,  were  applied  in  order  to  obtain  high  relia¬ 
bility^)  estimates  at  a  specified  lifetime.  In  the  analysis  the  unreliability  R  is  initially  determined 
from  application  of  a  discrete  joint  probability  density(JPD)  function.  The  function  provides  proba¬ 
bilities  associated  with  the  simultaneous  occurence  of  both  a  spectrum  load  and  a  material  strength 
value.  In  the  analysis  both  the  spectrum  load  scaling  factor  a  and  the  material  fatigue  strength  s, 
where  the  mean  fi,  =  Se,  are  represented  by  a  normal  PDFs  as  shown  in  Figure  1.  This  represen¬ 
tation  allows  for  application  of  the  JPD  analysis  in  addition  to  providing  for  potential  variability 
in  loading  and  material  strength.  The  scaled  version  of  the  spectrum  load  (o5t)  is  only  involved 
in  computing  N  in  Equation  1  but  not  the  R  estimate.  The  R  compututation  involving  only  a 
and  s  in  the  JPD  computation  is  therefore  simplified.  Substitution  of  a  for  oSl  is  valid  since  both 
are  normally  distributed  and  their  probability  computuations  are  independent  of  location.  That  is, 
the  scaled  and  unsealed  version  of  the  load  will  provide  identical  probability  estimates  with  respect 
to  the  JPD  computation.  Since  the  “event”,  of  identifying  a  particular  value  of  a  with  a  compo¬ 
nent,  is  independent  of  the  “event”  of  using  a  component  with  particular  fatigue  strength,  the  joint 
probability  that  a  =  a,  and  s  =  s,  occurs  simultaneously  can  be  written  as, 


P{a  =  ati,s  =  Sj)  =  P(a  =  q<)  •  P(s  =  Sj),  (5) 

where  i  =  1,2,3,  ...,«i  and  j  =  1, 2, 3, . . . ,  n2.  The  «i  and  n2  represent  the  number  of  events 
for  load  and  strength  respectively.  The  regions  Ac  and  As  where  the  events  occur  which  produce 
higher  probability  of  failure  are  bounded  hv  the  normal  PDFs  as  shown  in  Figure  1.  The  load  and 
strength  functions  are, 


(6) 


and 


/.(*)  = 


,v/2T 


a  2  V  a  *  ) 


(7) 


where  (/ia,<ra)  and  (/t5,<74)  are  the  population  means  and  standard  deviations  for  the  load  and 
strength,  respectively.  Referring  to  Figure  1  and  Equation  5,  the  JPD  can  be  written  as 


where, 


P  =  P  ■  P 

1  i,j  1  a,  ‘  s,i 


(8) 


Pa,  =  A  Qi  /a(Q,),  (9) 

PS]  =  Asrfa{s~),  (10) 

and  i  =  1,2, 3,  ...,ni  and  j  =  1,2,3, ...  ,/i2.  After  determining  the  joint  probability  values 
from  Equation  8,  the  corresponding  a  and  s  associated  with  these  probabilities  are  introduced  in 
Equations  2  and  1  respectively.  Thir  determines  a  specific  number  (A.j)  cycles  to  failure  of  the 
material  for  the  corresponding  probability  {Pij)  of  the  joint  occurance  of  a*  and  Sj.  The  lifetime 
estimate  Np  from  Equation  3  for  the  joint  o,  and  sj  event  is  obtained  from  the  following  application 
of  the  spectrum  load  data  {.$'/,( &)}  ^ ,  {Sm(k)}™'  and  { n(k )}fK  in  Table  1,  where  NK  is  the  number 
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of  spectrum  loads.  The  damage  fraction  for  a  specified  event  can  be  determined  from  Equation  4 
and  written  as 


n(k) 


iW 

- £  ",<*> 


(in 


where  n(k)s  are  the  spectrum  load  cycles  corresponding  to  the  original  tabulated  loads  Si(k).  From 
Equations  3  and  8,  the  lifetime  values  are  then  computed  from 


Np(ij)  =  1  /DFa 


(12) 


These  values  correspond  to  the  joint  probability  that  a  =  a*  and  s  =  Sj.  The  above  process  is 
repeated  M  =  n\n2  times, where  nx  and  n2  represent  the  number  of  mesh  points  associated  with  the 
tail  region  of  the  PDFs  in  Figure  1.  All  combinations  of  i  and  j  are  introduced  in  order  to  obtain 
paired  PXj  and  Np(ij)  values.  Ordering  only  the  Np  values  from  the  smallest  to  largest  and  retaining 
their  original  corresponding  P,j  probabilities  describes  a  discrete  PDF  representing  the  component 
probability  of  failure  Pj(t)  as  an  array  of  lifetimes  {Np(t)}^ ,  where  t  is  an  integer  defining  the 
ordering  of  the  Np  values.  See  a  graphical  display  of  a  PDF  in  Figure  2.  In  order  to  obtain  the 
unreliability  R  for  a  given  tx  in  Figure  2,  a  cummulative  density  computation  is  required.  This  is 
accomplished  by  selecting  an  ordered  value  from  Np  and  computing  the  sum 

=  (i3) 

i=l 

Note,  the  reliability  R  can  be  obtained  from  R  =  1-  R  and  the  lifetime  values  can  be  determined 
from  a  given  R. 


CONTAMINATED  PROBABILITY  DENSITY  FUNCTIONS 


In  order  to  illustrate  the  sensitivity  of  high  reliability  calculations  to  small  deviations  from  the 
assumed  models,  the  approach  taken  in  Reference  5,  is  applied.  Consider  the  situation  where  with 
a  high  probability  of  1  —  e,  samples  are  obtained  from  a  primary  PDF,  while  with  probabilty  e 
samples  come  from  a  secondary  PDF.  This  bimodal  probability  model  is  a  type  of  a  general  class 
referred  to  as  a  contaminated  models.  The  secondary  component  is  called  the  contamination  and 
the  probability  c  is  the  amount  of  contamination.  An  example  may  help  clarify  this  idea.  Consider 
the  situation  where  99%  of  the  specimens  are  obtained  from  a  population  of  “good”  specimens  while 
the  remaining  1%  of  the  time  consistently  lower  strength  measurements  are  obtained,  either  due 
to  manufacturing  defects  or  to  faulty  testing.  The  primary  PDF  would  correspond  to  the  “good” 
specimens,  the  contamination  would  represent  the  distribution  of  flawed  specimens,  and  the  amount 
of  contamination  is  e  =  0.01.  The  following  procedure  is  introduced  in  order  to  examine  the  effects 
of  computing  high  reliability  values  when  uncertainties  exist  in  selecting  the  PDFs  for  the  joint 
density  computations.  Initially,  values  are  obtained  from  the  JPD  computation  using  PDFs  fa  and 
/,  in  Equations  6  and  7.  Another  R  value  is  then  obtained  by  applying  the  PDFs  with  a  small 
amount  of  contamination  e. 

The  fa  PDF  with  variance  contamination  for  the  strength  data  is  written  as, 

fet  =  (1  -  < )/.(/«.,  a])  +  «/»(/<*,  *i2o?)i  (14) 
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where  pts  and  <7,  are  defined  in  Equation  7,  kx  is  a  scaling  factor  and  100«  is  the  percent  contami¬ 
nation.  A  similar  contaminated  distribution  for  fa  representing  load  a  can  be  written  as 

fee  =  (1  -  (15) 

with  scaling  factor  Aq  and  /ia  and  <ra  obtained  from  Equation  6.  Variance  contamination  produces 
effects  which  can  be  considered  to  represent  uncertainties  associated  with  the  selection  of  competing 
PDF  models. 

A  strength  distribution  with  mean  location  contamination  is 

ft,  =  (1  “  0/.(/*.,*J)  +  «/•(/*•  ±  *2<r„<V2),  (16) 

where  k2  is  a  scaling  factor  for  p,  and  the  sign  determines  which  tail  of  the  function  is  to  be 
contaminated  and  a/2  is  the  variance  for  /i,  ±  k2as.  The  contaminated  function  for  load  a  can  be 
determine  in  a  similar  manner.  The  location  contaminated  PDF  can  represent  the  rare  occurance 
of  exceptionally  high  loads  or  the  unusually  low  material  strength  of  a  degraded  or  non-conforming 
component  in  computing  the  reliability.  For  t  =  .01  and  k\  =  4,  graphical  results  in  Figure  3,  show 
an  almost  undetectable  difference  between  the  original  normal  PDF  and  the  contaminated  one.  A 
linear  relationship  to  obtain  R  from  the  JPD  function  application  can  be  obtained  by  combining 
both  contaminated  and  uncontaminated  functions  such  that, 


R*  —  (1  ~  ci)(l  —  £2)^00  +  Ml  —  f2)^io  +  £2(1  —  ei)Roi  +  (17) 

where  lOOej  and  100e2  are  the  percent  contamination  in  the  a  and  s  distributions  and  Rmn  represents 
reliabilities  obtained  from  contaminated  conditions  designated  by  m  and  n.  If  m,n  =  0,  then  there 
is  no  contamination.  If  m  =  1  and  n  =  0  then  fa  is  contaminated.  If  m,n  =  1,  then  both  fa  and  /, 
are  contaminated.  For  example,  if  there  is  contamination  of  the  strength  PDF  with  respect  to  the 
variance  then, 

/?"  =  ( 1  —  e2)/?oo  +  t2^oi-  (18) 

The  Rmn  values  are  obtained  from  R  in  Equation  13.  This  procedure  provides  an  effective  approach 
for  demonstrating  the  effects  of  PDF  uncertainties  in  determining  high  R  values. 


MODIFYING  TAILS  OF  THE  PDFS 


A  modification  of  the  fa  and  f,  PDFs’  upper  and  lower  tail  regions  respectively  was  introduced 
in  the  analysis  to  investigate  truncation  effects  as  suggested  in  Reference  3.  A  proposed  modifica¬ 
tion7  is  shown  in  Figure  4.  The  lumping  method  of  truncation  shown  in  the  figure  was  selected  so 
that  the  area  under  the  modified  PDF  remains  equal  to  one.  This  was  accomplished  by  determining 
the  area  under  extreme  tail  regions  associated  with  the  probabilities  P0ni  and  P,n7  obtained  from 
Equations  20  and  21.  These  areas  were  lumped  at  the  truncation  points  zt  and  z2  for  a  and  s  and 
the  reliability  Ri  values  were  determined  from  Equation  13,  with  the  lumped  /0,  /,.  A  comparison 
was  then  made  between  lumped  and  unlumped  results  in  order  to  determine  if  the  effects  of  the 
uncertainties  in  the  extreme  tails  are  significant  in  computing  high  R  values.  The  modification  also 
represents  a  substantial  difference  in  the  lower  tail  region  when  compared  with  the  original  PDFs. 
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The  R  computation  using  the  lumped  PDFs,  in  Figure  4,  is  the  same  as  that  described  previously, 
except  in  Equation  S, 


where, 


p  -  p  .  p 

1  n\nq  —  *  on,  1  s„2  i 

(19) 

too 

Pcni  =  fa-  da, 

(20) 

Jzi 

[*7 

Psn2  =  /  fs  ■  ds , 

J—OO 

(21) 

and  Z\  =  na  +  AjCTc  z2  =  n,  +  A2<r,,  are  the  truncation  points  of  the  PDFs  shown  in  Figure  4.  The 
fa  and  /,  PDFs  are  defined  in  Equations  6  and  7.  The  R  computation  procedures  are  then  applied 
using  the  newly  defined  PQni  and  PSnj  values. 


PDF  PROBABILITY  AT  TRUNCATION  POINTS 


The  lumping  procedure  described  previously  was  introduced  in  order  to  determine  the  effects 
on  the  R  values  from  modifications  to  the  extreme  tail  of  the  PDFs. 

In  applying  Equations  20  and  21,  it  is  assumed  that  the  PDFs  and  the  truncation  points  Z\ 
and  z2  are  known  exactly;  which  is  usually  not  the  case  for  engineering  problems  involving  material 
strength  or  loading  measurements.  Since  the  accuracy  in  estimating  the  truncation  point  locations 
is  essential  in  determining  the  importance  of  correct  extreme  tail  representation  m  obtaining  high 
R  values,  the  following  study  was  performed  involving  determination  of  the  reliability  and  quantile 
values  at  selected  truncation  points  as  a  function  of  sample  size.  Quantile  values  of  /,  can  written 
as, 

S„  =  /*,  -  K  •  or,,  (22) 

where  and  at  are  known  mean  and  standard  deviation  of  /,.  The  K  values  can  be  selected  to  define 
points  where  truncation  may  be  introduced  in  the  high  reliability  computation.  Unfortunately,  the 
and  <r,  values  are  not  known  sufficiently  well  for  an  accurate  measurement  of  Sq  unless  very  large 
data  sets  are  applied. 

The  following  simulation  process  examines  the  sample  size(n2)  effects  in  computing  the  trun¬ 
cation  points.  In  the  simulation  process,  a  »2  set  of  s(  normally  distributed  values  are  selected 
from 

*i(0  =  /<*0  +  ^  ■  Q,),i  =  1,2,3,...  ,n2  (23) 

where  the  Qi  values  are  obtained  from  a  standard  normal  distribution  with  p,  =  100  and  v,  =  .10. 
The  v,  value  is  the  CV  =  and  /<,  is  the  population  mean.  From  the  .«t  values  the  mean  st 

and  variance  VAR(sj)  are  determined.  An  estimate  of  the  population  quantile  Sq  is  then 

Sq  =  st-K.(VAR(st))1'2  (24) 

The  probability  Pj,  where  Pf  equals  Prob[.st  >  5V]  can  be  estimated  by  the  proportion  of  st  values 
which  are  greater  than  Sq.  The  process  involving  Equations  23  and  24  is  repeated  many  times  so  that 
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the  Pj  and  Sq  values  are  not  effected  by  further  increasing  the  number  of  simulations.  The  range  of 
Pt  and  Sq  values  is  a  measure  of  the  statistical  stability  of  the  sampling  process.  Particular  quantile 
values  such  as  1%  or  99%  can  be  obtained  by  forming  a  cummulative  probability  distribution  for  Sq 
and  Pj.  The  probabilities  are  ordered  from  the  smallest  to  largest  and  their  percent  is  determined 
from  their  numbered  position  in  the  ordering  which  is  divided  by  the  total  number  of  simulations. 
The  Sq  quantile  is  the  value  at  the  same  numbered  position  of  interest. 


RESULTS  AND  DISCUSSION 


In  this  section,  results  are  obtained  for  the  AIIS  round  robin  problem  in  Reference  2.  A  thin 
AISI  4340  steel  plate  with  a  central  hole  is  loaded  by  the  Felix  28  spectrum.  The  ASTD  S  -  N  curve 
coefficients  are  used:  C  =  3.5x10*’,  D  =  -1.47164  and  Sp  =  54.5KSI.  The  Su  value  is  180KSI  in 
Equation  2.  The  number  or  mesh  points  and  n2  of  the  PDFs  in  the  JDF  computation  (Equation 
5)  are  each  50,  where  Ca  =  /<0  and  C,  =  /r3( see  Figure  1). 

In  Figure  5,  representation  of  reliability  as  a  function  of  life(Ap)  is  shown  for  selected  CV 
values  used  in  defining  the  PDFs  fa  and  f,  in  the  R  computation,  with  a  mean  load  factor(/rQ) 
of  .70.  The  results  show  a  very  rapid  initial  decrease  in  reliability  followed  by  a  more  moderate 
decline  in  reliability  as  the  lifetime  values  increase.  Results  were  the  same  for  other  fia  values.  R 
estimates  also  decreased  with  an  increase  in  the  assumed  CV  value.  For  example,  a  CV  =  .05  for 
both  PDFs  provides  a  much  greater  R(.9(i,))  estimate  than  for  a  CV  =  .07  which  is  .9(6)  when 
Np  =  100.  The  results  designated  by  the  *  were  obtained  from  applying  lumped  PDFs  for  A,  = 
\2  =  3.5.  These  values  of  A  are  almost  the  maximum  amount  of  truncations  that  will  provide  the 
R  =  .96  requirement.  The  maximum  truncation  was  avoided  because  the  R  =  .96  would  be  met  at 
the  discrete  probability  value  of  the  discrete  joint  distribution  which  includes  the  lumped  values. 
The  R  =  .96  value  from  the  lumped  PDF  differs  slightly  from  the  unlumped  PDF  where  R  =  .9,<;4 
when  Np  =  50  and  the  CV  =  .8.  This  result  indicates  that  modifications  of  the  extreme  lower  tails 
of  the  PDFs  do  not  cause  large  differences  in  computing  high  reliability. 


Figure  6  shows  reliability  as  a  function  of  life  for  selected  fia  load  factors.  As  in  Figure  5,  there 
is  a  reduction  in  reliability  with  an  increase  in  Np  values.  There  is  also  an  obvious  decrease  in 
reliability  with  an  increase  in  load  factor.  In  the  case  where  Np  =  275,  an  /r0  of  .5  and  .6  resulted 
in  R  =  .9n  and  .96  respectively.  This  shows  a  substantial  sensitivity  in  computing  R  for  a  relatively 
small  differences  which  may  occur  in  estimating  /ia. 

Figure  7  also  shows  the  reliability  values  as  a  function  of  life(Ap).  The  assumed  CV  value  for 
the  PDFs  is  .07  and  applied  loads  are  fia  =  .5  and  .7.  The  dash  lines  represent  results  from  applying 
the  contaminated  PDFs  described  in  Equations  14  and  15,  for  the  case  of  1%  contamination  in  both 
PDFs  and  k i  =  4.  This  almost  undetectable  level  of  contamination  caused  a  drastic  reduction  in 
reliability  for  R  >  .9C.  The  results  from  the  contaminated  PDFs  application  which  represent  the 
uncertainties  in  assuming  a  specific  function,  demonstrate  the  importance  of  identifying  precise 
PDFs  in  computing  component  high  statistical  reliability.  In  contrast,  at  Np  =  130,  the  unlumped 
result  for  of  .952  is  only  moderately  reduced  from  the  value  of  .96  obtained  for  lumped  PDFs  for  Ai 


and  A2  =  3.5. 


In  Table  2,  reliability  results  are  tabulated  using  both  lumpedf/?^)  and  unlumped(/?t/i)  PDF 
applications  for  selected  CV  values.  These  results  show  a  reduction  in  Rl  and  Rul  values  with 
increasing  CV  values.  In  the  case  of  CV  =  .05  and  .06,  the  Ri  estimate  of  .976  was  the  maximum 
obtainable  because  of  the  discrete  nature  of  the  PDF  truncation  procedures.  Comparing  the  values 
of  Ri  —  .976  with  Rul  =  -96  for  CV  =  .07  shows  a  relatively  small  difference  in  the  reliability 
values.  This  result  indicates  that  uncertainties  in  modeling  the  extreme  lower  and  upper  tails  of  the 
PDFs  cause  relatively  small  differences  in  computing  high  R  values.  The  issue  that  is  important 
involves  the  substantial  reduction  in  the  Ri  and  Rul  values  with  increasing  CVs.  Since  the  CV 
values  are  often  estimated  from  coupon  data  that  are  assumed  to  be  relevent  to  actual  component 
behavior,  substantial  uncertainties  can  result  in  estimating  R. 

Table  3  shows  effects  in  computing  R  from  applying  both  contaminated  and  uncontaminated 
PDF  applications.  In  the  case  where  fia  =  .5  and  Np  =  80,  the  uncontaminated  PDF  result  is 
R  -9i2.  When  there  is  a  1%  contamination  of  PDF  for  the  load,  the  value  is  reduced  to  .96. 
Contaminating  the  strength  PDF  by  1%  resulted  in  a  reduction  from  twelve  nines(uncontaminated) 
to  three  nines(contaminated)  in  the  R  estimate.  Very  similar  results  wereobtained  for  contamination 
of  both  PDFs.  As  previously  shown  in  Figure  7,  obtaining  extremely  high  reliability  greater  than 
twelve  nines  will  not  provide  the  necessary  conservative  estimate  for  R  if  there  is  even  a  very  small 
amount  of  uncertainty  in  assuming  a  PDF  in  the  R  computation.  The  case  where  =  .7*  represents 
application  of  both  the  lumped  and  contaminated  PDFs  in  the  R  computation.  The  lumped  PDF 
result  without  contamination  showed  R  =  .9T6  which  is  reduced  to  ,928955  when  the  contaminated 
PDF  was  also  applied.  Again,  this  substantial  reduction  in  R  demonstrates  the  importance  in  the 
accuracy  of  the  PDF.  Potential  uncertainties  associated  with  defining  the  extreme  tails  of  the  PDFs 
become  insignificant  relative  to  the  accuracy  of  PDF  assumption  in  computing  high  R  values.  The 
table  also  shows  that  by  increasing  the  CV  value  the  Np  value  is  reduced  but  the  reduction  in  R 
from  the  contaminated  PDF  are  the  same  as  those  for  CV  =  .07.  Summarizing  the  results  in  the 
table:  it  is  critical  in  computing  high  statistical  R  values  that  the  PDFs  are  known  almost  exactly 
while  uncertainties  in  the  extreme  tails  of  the  PDFs  are  relatively  insignificant. 

Table  4  provides  results  similar  to  t  hose  in  Table  2  except  that  the  load  p0  is  varied  in  order 
to  examine  the  effects  of  uncertainty  in  determining  the  mean  scaling  load  factor  in  the  reliability 
-  life  estimating  process.  The  results  again  show  a  substantial  difference  in  R  for  an  uncertainty  in 
fia.  For  example,  when  /ia  =  .7,  R  =  .9(i  and  for  fiQ  =  .8.  R  =  .93662  which  implies  that  there  will 
be  one  failure  and  338  failures  in  a  million  for  loads  fi„  =  .7  and  .8  respectively.  This  substantial 
difference  relative  to  uncertainties  in  estimating  indicates  that  the  R  computations  are  very 
sensitive  to  uncertainties  in  the  load.  The  table  shows  little  difference  between  Rul  and  Rl  for 
Po  >  -7.  When  /t0  <  .7  no  quantitative  comparisons  are  possible  because  of  the  PDF  truncations. 
The  R  value  for  p0  =  .7  shows  a  substantial  decrease  in  Rul  from  approximately  .96  to  .9289  for 
the  uncontaminated  and  contaminated  PDFs  respectively.  A  similar  result  is  shown  for  Rl  at  /iQ  = 
.7.  The  Rul  and  Rl  results  at  /<rt  =  .7  showed  a  small  effect  on  R  with  a  substantial  modification 
of  the  extreme  tails  of  PDFs.  This  shows  that  the  PDF  assumption  is  critical  in  computing  R  while 
accuracy  in  representing  the  extreme  tail  of  the  PDFs  is  much  less  critical. 
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In  Table  5  are  the  results  of  sampling  a  generic  normal  PDF  to  examine  the  stability  of  the 
statistical  estimates  at  the  potential  truncation  points.  The  median  probability  and  quantile  values 
are  shown  for  a  range  of  sample  sizes^i^).  Included  in  the  results  are  the  upper  and  lower  bounds 
on  the  98%  confidence  interval  on  the  median  estimates.  Reliability^^)  values  are  obtained  at  K  = 
3.5  and  4.75  where  truncation  may  introduced  in  high  reliability  computations  in  studying  effects  in 
PDF  tail  modifications.  This  was  done  in  order  to  examine  if  there  is  instability  in  Pj  at  the  points 
due  to  the  sample  size.  The  sampling  trials  were  repeated  6000  times  which  was  sufficient  to  ensure 
that  the  tabulated  values  would  not  change  with  additional  trials.  Results  from  the  table  indicate 
that  relatively  unstable  Pt  values  will  be  obtained  for  even  a  sample  size  of  100.  In  this  case,  the 
true  Pj  value  is  .93767  for  I\  =  3.5  but  the  simulation  result  shows  an  inner  confidence  range  of 
.928172  to  .9482  for  Pj  values  associated  with  5,  estimates.  This  uncertainty  in  the  Sq  location, 
reduces  the  validity  in  assuming  that  if  lumping  a  PDF  does  not  cause  a  substantial  change  in  R, 
then  the  PDF  will  be  adequate  for  computing  high  R  values  such  as  .96.  Another  more  obvious 
example  is  the  case  where  n 2  =  6  which  shows  an  intervals  of  .835123  to  .998  for  K  =  3.5  and 
.918483  to  1.0  when  K  =  4.75.  The  inner  9S%  ranges  of  5,  quantile  values  for  K  =  4.75  and  K 
=  3.5  also  show  a  substantial  overlap.  This  case  shows  that  even  if  the  lumping  process  provides 
results  showing  small  differences  in  R  between  K  =  3.5  and  4.75  (using  an  unverified  assumed  known 
truncation  point),  the  inference  is  meaningless.  That  is,  the  substantial  uncertainty  associated  with 
computing  R  at  unverified  truncation  points  prevents  making  any  assessment  regarding  the  need 
for  accurate  representation  of  the  extreme  tail  of  the  PDF  in  computing  .96  R  values. 

These  results  are  consistent  with  results  of  truncations  of  normal  PDFs8  where,  for  truncations 
of  less  than  10  to  20  percent  of  the  population,  quantiles  would  fall  within  permissible  limits  of 
random  variation,  unless  sample  sizes  are  very  large.  Reference  5  shows  various  levels  of  uncertainty 
associated  with  computing  high  reliability  from  a  stress-strength  statistical  model  as  a  function  of 
sample  size.  These  results  relate  directly  to  the  sample  size  issue  discussed  in  this  report. 

The  substantial  sensitivities  of  R  in  each  of  the  figures  and  tables  relate  to  uncertainties  in  only 
one  parameter,  while  the  others  are  held  constant.  In  design,  the  uncertainties  in  more  than  one 
parameter  such  as  fia  and  CV  could  cause  increased  R  sensitivity.  There  are  many  complex  issues 
involved  in  obtaining  a  component  population  PDF  for  effective  load  severity  scaling  parameter, 
over  lifetime.  There  is  no  industry  standard  approach  to  characterizing  the  load  history  and  limited 
experience  in  determining  loading  PDFs.  Therefore,  the  substantial  influence  on  high  reliability 
caused  by  loading  PDF  uncertainties  could  cause  a  serious  problem  in  the  implementation  of  a  high 
reliability  requirement. 

These  results  are  based  on  a  single  S-N  curve,  in  contrast  to  the  AIIS  round  robin  problem 
in  which  each  participant  used  a  different  .S-N  curve.  Thus,  very  substantial  variability  in  the  R- 
lifetime  relationship  can  be  expected  even  when  the  S-N  curve  shape  and  mean  fatigue  limit  stress 
is  fixed. 

These  results  and  the  previous  analyses  of  contaminated  PDFs  in  Reference  5,  support  the 
concern  expressed  in  Reference  4,  for  the  issue  of  a  decrease  in  reliability  caused  by  degraded  or  non- 
conforming  components.  The  approach  of  attempting  to  obtain  statistically  very  high  R  values(.9i2) 
to  compensate  for  uncertainties  in  assuming  a  PDF  may  not  provide  an  effective  margin  of  safety 
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or  conservatism  relative  to  a  .1)G  requirement. 


The  comparisons  between  Ri  and  Rui  do  not  directly  relate  to  PDF  modeling  for  design.  In  the 
approach  used  in  this  report  the  lumped  value  in  the  PDF  is  made  exactly  equal  to  the  extreme  tail 
of  a  known  PDF  with  which  it  is  being  compared.  In  the  design  process,  the  difference  of  interest  is 
between  a  truncated  assumed  PDF  and  a  “correct”  PDF  which  is  unknown.  The  lumping  approach 
used  in  this  study  would  tend  to  minimize  the  difference  between  Rl  and  Rm  relative  to  an  actual 
design  process. 

No  conclusion  can  be  reached  about  an  acceptable  degree  of  truncation  from  this  study.  For 
\iand\2  <  3.5  it  appears  that  truncation  is  not  acceptable  for  the  .96  requirement.  Variation  in 
R  from  less  than  one  “nine”  to  values  approaching  two  “nines”  were  obtained  for  idealized  condi¬ 
tions  which  minimize  reliability  differences  as  noted  previously.  For  X^andX^  >  3.5,  the  sampling 
results  indicate  that  it  does  not  appear  to  be  feasible  to  obtain  satisfactory  representation  of  PDF 
unless  very  large  data  sets  are  available.  More  important,  the  issue  of  acceptable  degree  of  trun¬ 
cation  appears  to  be  of  secondary  importance  relative  to  the  sensitivity  of  high  R  to  the  expected 
uncertainties  in  assuming  a  specific  PDF  representations. 


CONCLUSIONS 


Unstable  high  statistical  reliability  values  for  a  fatigue  loaded  component  can  result  from  un¬ 
certainties  in  assuming  the  PDF  model  and  determining  its  parameters  without  using  very  large 
data  sets  in  the  analysis. 

Estimates  of  the  extreme  tail  quantiles  and  their  corresponding  reliabilities  can  be  unstable 
unless  large  data  sets  are  used. 

Analysis  of  the  effects  of  extreme  tail  modification  does  not  provide  decisive  information  on  the 
adequacy  of  PDF  modeling.  Tail  modification  effects  on  reliability  are  small  relative  to  the  effects 
of  uncertainties  in  assuming  a  PDF  model. 

The  primai  i/  conclusion,  from  the  analytic  evaluation  in  this  report ,  is  that  computation  of  high 
statistical  reliability  may  have  little  or  no  association  with  actual  component  reliability. 
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Table  1:  Rainflow  low-high  load  sequence  derived  from  Felix-28 


k 

k  Sl  Sm  n(k) 

1 

2.80 

25.59 

354 

26 

42.63 

29.21 

207 

2 

2.80 

32.83 

334 

27 

42.63 

36.45 

1274 

3 

6.42 

29.21 

416 

28 

46.25 

21.97 

274 

4 

10.04 

29.21 

609 

29 

46.25 

25.59 

6239 

5 

10.04 

36.45 

1228 

30 

46.25 

29.21 

4274 

6 

10.04 

40.07 

810 

31 

46.25 

40.07 

604 

7 

13.66 

36.45 

2 

32 

49.87 

3.86 

268 

8 

17.28 

18.35 

140 

33 

49.87 

25.59 

956 

9 

17.28 

32.83 

78 

34 

49.87 

29.21 

2179 

10 

20.91 

32.83 

2061 

35 

53.49 

25.59 

2 

11 

20.91 

36.45 

90 

36 

53.49 

29.21 

116 

12 

24.53 

-7.00 

140 

37 

57.12 

25.59 

5 

13 

24.53 

18.35 

140 

38 

57.12 

29.21 

185 

14 

24.53 

36.45 

2040 

39 

60.74 

29.21 

25 

15 

28.15 

29.21 

833 

40 

64.36 

25.59 

7 

16 

31.77 

25.59 

346 

41 

64.36 

29.21 

8 

17 

35.39 

25.59 

7904 

42 

64.36 

32.83 

75 

18 

35.39 

29.21 

56 

43 

67.98 

29.21 

9 

19 

35.39' 

32.83 

71072 

44 

71.60 

29.21 

16 

20 

39.39 

43.69 

2529 

45 

75.22 

25.59 

7 

21 

39.01 

21.97 

3014 

46 

78.84 

18.35 

5 

22 

39.01 

25.59 

42825 

47 

78.84 

25.59 

1 

23 

39.01 

29.21 

6393 

48 

82.46 

21.97 

128 

24 

39.01 

43.69 

252 

49 

82.46 

29.21 

16 

25 

42.63 

25.59 

480 

50 

89.70 

25.59 

8 

12 


Table  2:  Reliability  from  lumped  and  un¬ 
lumped  PDFs  as  a  function  of  the  CVs 


CVs 

Rl 

Rul 

>  ,976 

•9n 

ISslI 

>  .976 

•97 

0.07 

.976 

•96 

0.08 

•954 

.9483 

0.09 

.93885 

0.10 

.93646 

.93533 

Rl,  -  Reliability  from  lumping  PDFs 
Rul  -  Reliability  from  normal  PDFs 
CVs  -  Coefficient  of  variations 

Np  =  80  Passes,  fia  =  .70  Load  factor,  Se  =54.5KSI  Strength 


Table  3:  Effects  of  individual  PDF  contaminations  on  reliability  estimates 


Ha 

!  Ru 

|  Rcl  | 

Res  | 

|  Rcls  | 

.5 

.07 

>  -9i2 

•96 

.93816 

.93814 

.7 

.07 

eIII 

>  .9e 

.93730 

.93186 

.928914 

a 

.07 

80.0 

>  .976 

.93733 

.93224 

.928955 

a 

.10 

4.0 

>  -96 

.93848 

.93092 

.928940 

*  Both  tails  of  PDFs  are  lumped  at  3.5  <r 

Ru  -  Reliability  uncontaminated  PDFs 

Rcl  -  Reliability  contaminated(l%)  load(a)  PDF 

Res  -  Reliability  contaminated(l%)  strength's)  PDF 

Rcls  -  Reliability  from  both  load  and  strength  PDFs  contamination 


13 


Table  4:  Reliability  from  lumped  and  un¬ 
lumped  PDFs  VS.  Load  fta 


f^atpka 

RuL 

Rl 

1ES9I 

>  -9l2 

>  .9r 

>  .9i2(.999814j* 

.97130 

■iW 

.9io 

.9r297 

0.7 

.9g[.998914]’ 

.9t[.998955]* 

0.8 

.999662 

.999724 

.984984 

.984710 

1.0 

.855658 

.851946 

Rl  -  Reliability  from  lumping  PDFs 
Rul  -  Reliability  from  unlumped  PDFs 
*  Results  from  contaminated  PDFs 
Np  —  80  Passes  And  CVs  =  .07 


Table  5:  Confidence  interval  on  Pr  and  SQ  quantile 


_ 1 

Normal  PDF  ft  = 

=  100,  CV 

=  .10 

H 

Median 

K  =  3.5 
Lower  Bound 

Upper  Bound 

Median 

K  =  4.75 
Lower  Bound 

Upper  Bound 

6 

.999480* 

.835123 

.9978 

.9S58 

.918483 

1.00 

(90.82) 

(37.23) 

(55.10) 

(86.67) 

(16.56) 

10 

.999669 

.925734 

.9786 

•9581 

-9l3 

(66.00) 

(85.04) 

(43.75) 

(78.07) 

(24.78) 

20 

.999721 

.981327 

.9S84 

.998403 

•9io7 

(65.48) 

(78.92) 

(53.36) 

(70.56) 

(34.27) 

50 

.999745 

.995914 

■MHKSJ 

.9S9 

.999845 

■9881 

(65.36) 

(73.92) 

(52.81) 

(63.95) 

(41.14) 

100 

.998172 

.9482 

.9S9 

(65.13) 

(71.05) 

(59.01) 

(52.68) 

(60.52) 

1000 

.999765 

.999524 

.999890 

•9«7 

_ 

(65.02) 

(66.96) 

(63.07) 

imkhqi 

(49.99) 

•Pr  Probability 

**  Corresponding  5,  quantile  value 
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Solid  Line:  Uncontaminated  PDF 
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Figure  3.  Containinated(l%)  and  uucontaininated  normal  PDFs 


Load  &  Strength  CV  =  .07 
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Figure  7.  Reliability  vs.  Ni>:  Contaminated- Uncontaminatod  PDFs 


Load  fia  Values 
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Figure  6.  Reliability  vs.  Np  for  specified  fia  values 
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